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Abstract 

We present an exploratory lattice study oi B ^ it semileptonic decay form 
factors using the nonrelativistic lattice QCD for heavy quark with Wilson 
light quark on a 16^ x 32 quenched lattice at (3 = 5.8. The matrix elements 
are calculated at eight values of heavy quark mass in a range of 1.5 — 8 GeV 
and with three values of light quark mass. The l/m^ corrections to the ma- 
trix elements are found to be fairly small except for the spatial component 
proportional to the B meson momentum. We find that the dependence 
of the form factor f~^{(f') near q^^^ becomes much stronger for larger heavy 
quark mass, which may suggest the increase of a pole contribution. We per- 
form a model independent fit of the form factors and study whether the 
dependence is consistent with the pole contribution. Although the soft pion 
theorem predicts f^{qf^^^,J = fs/ fn in the chiral limit, we observe a significant 
violation of this relation. 
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I. INTRODUCTION 



The exclusive semileptonic decay B 7r(p)/z/ would become an important process to de- 
termine the Cabibbo-Kobayashi-Maskawa (CKM) matrix element \Vub\, when high statistics 
experimental data becomes available by future B Factories, since the lattice QCD simula- 
tion enables us to compute the relevant form factors from the first principle. There was, 
however, a difficulty in treating a heavy quark with mass mq on the lattice, because pos- 
sible systematic errors of order amg could become unacceptably large for a typical lattice 
spacing a accessible in present simulations. Thus the previous lattice calculations of the B 
meson semileptonic decay form factors involved an extrapolation in the heavy quark 
mass from the charm quark mass regime to the 6-quark mass assuming a heavy quark mass 
scaling law, which could introduce a potential systematic error. The nonrelativistic lattice 
QCD (NRQCD) is designed to remove such a large uncertainty based on a systematic 
l/mg expansion, and one can simulate the 6-quark directly at its mass value. In this paper 
we describe the lattice calculation of the B — > ttIu form factors using the NRQCD for heavy 
quark. We investigate the heavy quark mass dependence of the form factors, by taking the 
mass of the heavy quark to cover a range of 1.5-8 GeV. 

The hadronic matrix elements of the B ^ n semileptonic decay are expressed in terms 
of two form factors and /° as 

{7rmV,\B{p)) = (p + k - q^^^^^) r{q') + q,^^^^^nq% (1) 

where q^ = Pfj, — k^. The covariant normalization of the meson fields is employed in this 
paper: 

{M{p)\M{p)) = 2poi2nf6%p-jf). (2) 

From the condition that the matrix element is not singular at q^ = 0, the form factors satisfy 
f~^{0) = /°(0), and the kinematical end point q"^^^ = {rriB — mT^Y corresponds to the zero 
recoil limit, where the lattice simulation works most efficiently. 

This paper is organized as follows. In the next section, the NRQCD formulation and 
numerical method to calculate the matrix elements are summarized. We describe the details 
of our numerical simulation in Section where we point out the subtleties in extracting 
the form factors: the definition of the heavy meson energy, the choice of the two independent 
matrix elements to determine /"*" and /°, and the procedure of chiral extrapolation. We ex- 
plain what we think is the best procedure and study the uncertainty by comparing the result 
with those from other procedures. Physical implications of numerical results are discussed 
in Section |rV|. We study l/m^ dependence of the matrix elements and q^ dependence of 
the form factors. The prediction from the soft pion theorem is compared with our data. 
In Section |V|, we discuss the systematic uncertainties contained in this work. Section [V^ is 
devoted to our conclusion. 
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II. LATTICE NRQCD 



A. Lattice NRQCD action 

The lattice NRQCD has been extensively used for the investigations of the heavy-heavy 
systems and hadrons containing a single heavy quark ||^. It is designed to remove the 
large mass scale ttiq from the theory using the 1 / rriQ expansion and to reproduce the same 
results as of the relativistic QCD up to a given order of l/mq. In this work, we employ the 
lattice NRQCD action including the ©(l/mg) termsQ 



NRQCD 



X 



Q{x + 4)-{l-SH)Qix) 
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and Q{x) is the effective two component spinor field, which describes the heavy quark. A^^^^ 
denotes the three dimensional Laplacian, and B is the standard clover-leaf type chromo- 
magnetic field. 

This action generates the following evolution equation 



GQ(f,t=l) 

GnixA+l] 



1 Ho 

2n " 
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^HoJ GQ{x,t = 0), 
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for which we apply tadpole improvement procedure 



U^{x)/uo with uo 



(Trt/piaq/S)^/"^ p. To avoid the singular behavior of high frequency modes in the evolu- 
tion equation, the stabilizing parameter n is chosen to satisfy the condition |1 — Ho/2n\ < 1, 
which leads to n > 3/2mQ. From a viewpoint of the perturbation theory, further constraint, 
Ho/2n < 1, is necessary to avoid singularities in some of vertices derived from the action 
(§). This point is again discussed in the last part of this section in connection with our 
choice of n in the simulation and the perturbative calculation. 

The four component spinor field ip{x) of the relativistic QCD is expressed in terms of 
two component spinor field Q{x) as 



^This action differs from that we used in our previous study of which is organized to remove 

0{ah.QCD /mo) error, at the cost of simulation speed. 
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where A is the symmetric lattice covariant derivative, ©(l/mg) correction appears in the 
lower component of ip{x), which affects the heavy-light current. 



B. Correlation functions 

We employ the standard simulation technique to calculate the hadronic matrix elements 
of the semileptonic decay. We calculate the three-point correlation functions 

= ^^e-^^-^/e-(^^-^>^^ (0|OB(f/, t/)V;t(x„ 4)0t (0, i,)|0), (9) 

Xf Xs 

where Ob and O^^ are the interpolating operators for B and vr mesons, respectively, and 
= q'j^ip is the heavy-light vector current. In this work, we use the Wilson quark to 
describe light quark q{x). Here we denote the heavy-light and the hght-light pseudoscalar 
mesons as B and vr, respectively, regardless of their mass parameters k and mg for simplicity. 
For t/ ^ 4 ^ ^ the correlation function Eq. (|^) becomes 

Cl^\p,k;tj,t,A) - ^^^e-^^'^'^^^f-^'^^^^^^^ (10) 

where Eb{p) and E.,^{k) denote the energy of B meson and pion, respectively. The exponent 
Eqqip) is not the total energy but the binding energy of the B meson, because the heavy 
quark mass uiq is subtracted in the NRQCD. We use the local interpolating operators for 
both of B and vr, and 

Zb{p) = (0|OB(0)|i?(p)), Z^{k) = (0|O^(0)|7r(^)) (11) 

are their matrix elements. 

In calculating Eq. @ we vary tj with fixed ti and tg in order to find out the region 
where the correlation functions are dominated by the ground state. The fixed ts is chosen so 
that the pion two-point correlation function is dominated by the ground state, as is shown 



in Section |ITl[ To obtain EqQ^p), ET,{k), Zb^p) and Z.„{k), we also calculate the two-point 



correlation functions with a finite momentum 



C^b\f. tf, U) = E e-^i'-^f {OB{xf)OUx^)) - exp{~E,Q{p){tf - U)), (12) 

Ci^\k;tf,U) =J2e'''-'^H0A^f)0tix,)) ^ exp(-E^(A;)(t; - t,)). (13) 



2E^{k) 
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Combining Eqs.(|T^), (0), and ([I3D, one can easily see that the matrix element is expressed 

as 

{Bmvi\m).... = pEB{v)^2E^{k) f ^ C^^\p,k;tfM) (14) 

Zb{p) Cy{k]ts,ti) 

for tf S> 4 ^ ^5 where Z = Z/\/2E. As expressed in Eq. (0), we use the two-point 
correlation function itself to cancel the exponentially decaying factor of pion, while use 
the values of EqQ obtained by fits to cancel the B meson's. One reason of this asymmetric 
procedure is that the pion two-point function is constructed from the light quark propagator 
with a point source at = 4, which is what we used to calculate the three-point function 
(1), and then we expect the statistical fluctuation mostly cancels between @ and (p!3D, while 
for the B meson exponential function, such a cancellation is not expected. In addition, as 
we mention in the next section, the two-point correlation function of B meson with the 
point source (0) requires larger time separation to reach the plateau than the three-point 
function (|^), for which the heavy quark source is effectively 'smeared' at ts- 



C. Perturbative corrections 

To relate the matrix element in the lattice theory to that in the continuum QCD, operator 
matching is required. We have calculated the perturbative renormalization factor Zy^ for 
the vector current at one-loop level using the lattice perturbation theory 

V-^ = Zv, F;^" = Zy^ q^.ij, (15) 

where q is the Wilson light quark and ip is defined in Eq. (H). Zy is the ratio of the on-shell 
S-matrix elements in the continuum theory with MS scheme and that in the lattice theory. 
In our definition, Zy contains the leading logarithmic term, alog(mQa), which comes from 
the continuum renormalization factor. 

In calculating Zy we use the massless Wilson quark and the external momenta are taken 
to be zero. We did not take into account the one-loop operator mixing with higher derivative 
operators, since there are already 0{a) errors at tree- level from the Wilson quark action. 
The one- loop coefficient is modified with the tadpole improvement For the mean link 
variable we use Uq = except for the light quark wave function renormalization, 

for which we use another possible definition, uq = 1/8/tc [l^. Their one- loop perturbative 



expressions are used to determine the perturbative coefficients of Zy^. 
The results for the one-loop coefficient Cy in 



l+9^Cy^ (16) 



are presented in Table | for several values of {mQ,n). These values contain the leading 
logarithmic contribution, log (mg a) /47r^. The values of Zy^ with two choices of the lattice 
coupling constant gy{n/a) = 2.19 and gy^l/a) = 3.80 are plotted as a function of l/mg 
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in Fig. H]. We observe that the spatial component of the vector current receives larger 
perturbative corrections than the temporal one. On the other hand, the l/mq dependence 
is rather stronger for Zy^ than for Zy.. 

When we discuss the l/ms dependence of the renormalized matrix elements in Section 
P^ , we multiply the leading logarithmic factor 

/ ^ ^ \ 2/11 
0{ms/m^ ) - [ayimr'^) ) ^^^^ 

to cancel the logarithmic divergence in the infinite heavy quark mass limit due to the anoma- 
lous dimension of the heavy-light current. 

The perturbative correction for the heavy quark self-energy is also calculated, and the 
B meson mass is given through the binding energy of the heavy-light meson EgQ{p = 0) as 

niB = Zmmg - Eo + EgQ{p = 0), (18) 

where the energy shift Eq and the mass renormalization Zm are obtained perturbatively 

Zm = l+ g^B, (19) 
Eo = g'A. (20) 

The tadpole improved coefficients A and B are also given in Table |. 

For a historical reason, the stabilizing parameter we have used does not always satisfy 
the condition n > S/ttlq, which is necessary to avoid divergent tree level vertices, while the 
simulation itself is stable with the condition n > 3/2mQ. We, therefore, quote the results at 
tree level in the later sections as our main results. We estimate the size of the renormalization 
effect with the one-loop coefficients obtained with the combinations of mg and n, for which 
n's are larger than those we have used in the simulation and the perturbation theory exists. 
Although this estimation is certainly incorrect, it gives some idea for the one- loop effect, 
especially because the n-dependence of the simulation results is observed to be very small 
(Section UTTdI) . 



III. SIMULATION DETAILS 

In this section, we describe the numerical simulation in detail apart from discussions 
on physical implications of the results, which will be discussed in the next section. After 
summarizing the simulation parameters, two-point correlation functions of vr and B mesons 
with finite momenta are discussed. We describe how to extract the matrix elements and the 
form factors from the three-point correlation functions. Finally, the chiral extrapolation of 
the matrix element is discussed. 
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A. Simulation parameters 



The numerical simulations are performed on a 16'^ x 32 lattice with 120 quenched gauge 
configurations generated with the standard plaquette gauge action at /?=5.8. Each configu- 
ration is separated by 2000 pseudo-heat-bath sweeps after 20000 sweeps for thermalization 
and fixed to the Coulomb gauge. The Wilson quark action is used for the light quark at 
three k values 0.1570, 0.1585 and 0.1600, which roughly lie in the range [ms,2ms], and the 
critical hopping parameter is fi;c=0. 16346(7). The boundary condition for the light quark 
is periodic and Dirichlet for spatial and temporal directions, respectively. The light quark 
field is normalized with the tadpole improved form 
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^1 — ?>k/Akc according to [0. The 
tadpole improvement is also applied for both the NRQCD action and the current oper- 
ator with the replacement of Uf^/u^ using the average value of a single plaquette 
Mo = (Trf/,,,,/3)i/4 = 0.867994(13). 

The lattice scale is determined from the p meson mass as a~^=1.71(6) GeV, although 
we expect a large 0(a) error for nip with the unimproved Wilson fermion. The results for 
the TT and the p meson masses and the pion decay constant are summarized in Table 0. 
The heavy quark mass mq and the stabilizing parameter n used in our simulation are 



(21) 



where mq = 2.6 and 0.9 roughly correspond to b- and c-quark masses, respectively. 

For mq=2.1 and 1.2 we performed two sets of simulations with different values of n, 
though the statistics is lower (=60) for {mq,n) = (2.1,2) and (1.2,3). Since the different 
choice of n introduces the different higher order terms in a in the evolution equation, the 
choice of n should not affect the physical results for sufficiently small a. The small depen- 
dence of the numerical results on n is also crucial for our estimation of the perturbative 
corrections. 

The spatial momentum of the B meson (p) and the pion (k) is taken up to |A;| < 
-\/3-27r/16, which corresponds to the maximum momentum of ~ 1.2 GeV in the physical unit. 
We measure the three-point correlation function at 20 different momentum configurations 
(p, k) as listed in Table ^TT[ The momentum configurations which are equivalent under the 
lattice rotational symmetry are averaged, and the number of such equivalent sets are also 
shown in Table |ITI| . 

The light quark propagator is solved with a local source at ti=A, which is commonly 
used for the two-point and three-point functions. The heavy-light vector current is placed 
at tg = 14, which is chosen so that the pion correlation function is dominated by the 
ground state signal. The position of the B meson interpolating operator is varied in a range 
= 23 — 28, where we observe a good plateau as shown later. 
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B. Light-light meson 



In order to obtain the form factors reliably, it is crucial to extract the ground state of the 
B meson and the pion involving finite momentum properly. In Fig. ^ we show the effective 
mass plot of pions with finite momentum at k = 0.1570 and 0.1600. The spatial momentum 
k = {k^, ky, kz) is understood with the unit of 27r/16. This notation will be used throughout 
this paper. Although higher momentum states are rather noisy, we can observe a plateau 
beyond t = 14. We fit the data with the single exponential function to obtain the energy 
E^^ik) shown by the horizontal solid lines in Fig. ^ 

Figure ^ shows the energy momentum dispersion relation of pion, where the solid lines 
represent the relation in the continuum = + k"^. We observe a small discrepancy 

between the above relation and the data, which indicates the discretization error. However 
the disagreement is about 1-1.5 standard deviation and only a few percent. 



C. Heavy-light meson 

To compute the B meson two-point correlation functions, we employ the smeared source 
for heavy quark as well as the local source, with the local sink for both cases. The smearing 
function for the heavy quark is obtained with a prior measurement of the wave function 
with the local source. In Fig. ^ we plot the effective mass for both the local-local and the 
smeared-local correlation functions at mg = 2.6 and k = 0.1570, 0.1600. The plateau is 
reached beyond t = 16 for the local-local, while the smeared-local exhibits clear plateau from 
even earlier time slices. We obtain the binding energy with a fit range [16, 24] for both types 
of the correlation functions and for all momenta, and the results are consistent in all cases. 
The binding energy averaged over the results fitted from the local and the smeared sources 



are listed in Table |I^ together with the values in the chiral limit. In Table |V|, we also 
listed the binding energy for the vector meson B* measured with the local-local correlation 
function, which are used in later discussions on the B* pole contribution to the form factors. 
It is also worth to note that the values of Eqq obtained with different stabilizing parameter 
n is consistent with each other within their statistical errors. 

The dispersion relation for the B meson takes the following nonrelativistic form 

E,q{p) = EggiO) + -^f + 0(l/m|), (22) 

where the kinetic mass m^in should agree with the rest mass rriB (0) in the continuum limit. 
Since we use the NRQCD action correct up to 0(l/mQ), including higher order terms in 
l/rriB in Eq. (|2^ ) does not make sense. In Fig. ||, Eqqij)) is shown as a function of "p^ at 
mg = 2.6. The solid lines represent the relation (|2^ ) with m^in = uib determined through 
the tree level relation rriB = mq + EqQ{Q), which reproduce the data quite well. With the 
one-loop correction (ITS) the agreement becomes even better as presented with the dashed 
lines in the figure. 
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D. Three-point function and matrix elements 

Figure 1^ is the effective mass plot of the three-point function at tuq = 2.6 and k, = 0.1570, 
0.1600. The horizontal axis represents the time slice on which the B meson interpolating 
operator is put, and the vertical axis corresponds to the binding energy of the B meson. 
The horizontal solid lines represent the binding energy EgQ{p) determined from the two- 
point correlation functions. The figures display that the three-point correlation functions 
are dominated by the ground states beyond t = 23, and there they give the consistent values 
for Eqqip) with ones extracted from the two-point functions. Therefore, in this region we 
can use Eq. (JT^) together with the results of the two-point correlation functions to extract 
the matrix elements. 

It is useful to define the quantity as 

V^K.vMWf))^, (23) 

pEB{p)\j2E^{k) 

because it is defined only through the residue of the two- and three-point correlation func- 
tions without the knowledge how one defines the meson energies. Since there are uncertain- 
ties in the light-light and heavy-light meson dispersion relations, it is better to deal with 
the quantity which is free from the ambiguity. Moreover, is the quantity which has the 
infinite mass limit in the heavy quark effective theory. When the perturbative correction is 
incorporated, Q{mB /m^s'^''^) given by Eq. (|17D is multiplied to V^. Therefore V"^ is suitable 
quantity to study the l/mg dependence. 

For the spatial components of V'^, we also define the scalar products 

Up{p, k) = , Uk{p, k) = ^ . 24 

In Table we list the values of V4, Up, and Uk for all momentum configurations {p, k) at 
rriQ = 2.6 and k = 0.1570. In this table, we also list the values of determined with the 
tree level mass relation ([T8|) for the B meson. 

We have investigated the n-dependence of at mq = 2.1 with n = 1 and 2 and at 
rriQ = 1.2 with n = 2 and 3, using the first 60 configurations on which {mq^n) = (2.1,2) 
and (1.2,3) data are measured^. For both of the heavy quark masses we observed small 
dependence on n, which is at most 1%, 8% and 2% for V4, Up and Uk respectively, and 



^ We note that n-dependence should be studied on the same configurations. In some of the figures, 
there appear large deviations for the data with difi'erent n but the same mq. However, in these 
graphs only the results for {mq,n) = (2.1,2), (1.2,3) are obtained from the first 60 configurations 
and the results for the other combinations of {mq,n) are obtained from the entire 120 configura- 
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smaller than their statistical error. In the present work, therefore, we regard them to be 
sufficiently small to estimate the size of the renormalization effect in the manner described 
in Section [11 C| . 



E. Form factors 

To convert V4, f/p, and Uk to the form factors, we need to assume certain dispersion 
relations for Eb{'P) and ET^{k). One method is to use the values obtained from the dispersion 
relation measured in the simulation. This, however, suffers from the large statistical error 
for the finite spatial momenta. Alternatively, we adopt the following relativistic dispersion 
relations for both the B meson and the pion. 

Eb{p) = ^ml+fP, E^{k) = ^ml + P, (25) 

where the measured rest mass is used for and m^. These relations are almost satisfied 
as shown in Figs. |^ and ^ for light-light and heavy-light mesons, respectively. 

Using the relations Eq. (|25|) , the form factors are easily constructed from V^. First, we 
calculate P{q^) with 

, , j2EB(p)\/2EJk) 

and f^{q^) is similarly obtained from (p -|- k)'^V^ substituting the value of /° determined 
above. 

For py^O and k 0, f'^ and are not uniquely determined from V4, Up, and Uk- In this 
case there is an additional relation among V^'s, which should be satisfied when the Lorentz 
symmetry is restored. For p k this relation reads 

EB{p)Up + E^{k)Uk = %. (27) 

We examine this condition for iq = 6,9, 14 and 16 {ig is referred in Table pn| ). Figure |^ 
compares LHS and RHS of Eq. ( pTf ) at k = 0.1570 for iq=Q, with the tree level dispersion 
relation for Eb- This figure exhibits a difference of about 15%. In other cases of iq, similar 
amount of the discrepancy is observed. The size of this systematic effect is consistent with 
the naive expectation for 0{a) error. 



tions. It seems that these large deviations seem to arise from the statistical fluctuation caused by 
the remaining 60 configurations for which there is no the data with (mQ,n) = (2.1, 2), (1.2, 3) . 
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F. Chiral extrapolation 



To obtain the form factors at the physical pion and B meson masses, it is necessary to 
extrapolate the results to the chiral limit. There is, however, still a subtlety in the chiral 
extrapolation, because the light quark mass dependence of the matrix elements or the form 
factors are not well understood. In principle, the chiral limit of the matrix elements or the 
form factors must be taken using the result of the chiral effective theory as a guide for its 
functional form. For the B tt semileptonic decay the heavy meson effective theory with 



chiral Lagrangian gives such an example |]Tl|-[lt 

At least the heavy meson effective theories tell us that the matrix elements or the form 
factors depend on v ■ k, where is the 4- velocity of the B meson. At the zero pion 
momentum, the quantity v ■ k could potentially give linear dependence in m^, which could 
result in a ^ym^ dependence. The zero recoil limit in the heavy meson effective theory gives 



the following relations for the matrix element and the form factor: 

{7rik = 0)\V,\Bip = 0)) = (mB + m^)fiqlJ=mB^, (28) 

Jtt 

Assuming the linear dependence of Jb, fn, and on rrig, at least in the zero recoil limit 
the matrix element should have linear dependence on m^. In the following analysis, we take 
the chiral limit of the matrix elements assuming the linear dependence on ruq in any case of 
(p, k), although there is no proof. 

Figure || shows the chiral extrapolation of the matrix element with the form 

{n{k)\V^\B{p)) =av + bvmg, (29) 

where rriq = 1/2k — 1/2kc. The data itself do not show any sign of nonlinear behavior at 
least around the strange quark mass. The form factors f^{q^) and P{q^) at the physical 
pion mass are extracted after extrapolating the matrix elements to the chiral limit using 
Eq.(Pi. 



IV. PHYSICAL IMPLICATIONS 

In this section we discuss the physical implications of our results, which include the I/ttib 
dependence of the B ^ tt matrix elements and the dependence of the form factors. The 
prediction from the soft pion theorem is compared with our data. 



A. l/uiB dependence 

The heavy quark effective theory predicts that the properly normalized 5 — > vr matrix 
element has a static limit, hence it can be described by an expansion in the inverse heavy 
meson mass I/tub whose leading order is a function of the heavy meson velocity = Pf^/niB, 
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(30) 



Similar arguments for the heavy-light decay constant suggested that the quantity /b^ttib 
has the static limit while numerical simulations have shown that the l/m^ correction is very 
large. On the other hand, the l/m^ dependence of the form factors have been studied only 
in the D meson region Therefore it is important to study the l/m^ dependence of 

the matrix elements at fixed values oi v ■ k. 

Except for p = 0, fixing p is not quite identical to fixing v ■ k, since the velocity changes 
depending on the heavy meson mass. Thus it is awkward to use the matrix elements with 
nonzero p. In the special case of p = 0, LHS of Eq. ( pO] ) is nothing but the matrix elements 
V4, Up and Uk, defined in Eqs. (p3D, and (p4D, multiplied by the rriB independent factor. 

In the following analysis, we confine ourselves to examine the following quantities for the 
sake of simplicity: 

/ (1) (2) \ 

V,{p = 0,k) = vP{l + ^ + ^ + ---], (31) 
y niB j 

( (1) (2) \ 

J7,(p = 0, fc) = ?7l°M 1 + ^ + \ + ■ ■ ■ , (32) 
Up{p = 0,k) =\im Up{p,k) (33) 

= — l + ^ + \ + , (34) 

rriB \ ruB m% j 

for which we explicitly show the form of the l/niB expansion. All of the coefficients in these 
expansions are a function of k. 

In Figs. H and |T0| we show the l/m^ dependence of V4 and f/^, respectively, at k = 0.1570. 



The l/rriB correction is not significant for these quantities and almost negligible around the 
B meson mass. This result exhibits a sharp contrast to the mass dependence of the heavy- 
light decay constant Jb^/'^b, for which the large l/rriB correction to the static limit is 



observed. Results of the linear and quadratic fit in l/m^ are listed in Table ^ for V4 and 
in Table [VT] for Uk. 

We note here that x^/dof are less than unity for most cases of V4, Uk, and also Up, which 
will be mentioned in the next paragraph, though they do not exactly judge the goodness of 
the fits for such data, which are correlated for different rriQ. 

In order to do the same discussion for Up, which is defined in the — > limit, we 



extrapolate the finite p results to the vanishing p point as shown in Fig. |TT[ There is little 
p^ dependence observed and we employ a linear extrapolation in p^. In Fig. |12| we plot mBUp 
as a function of l/m^ at k = 0.1570. In contrary to the other matrix elements we observe a 
sizable l/niB dependence. Table [VIII| summarizes the results of linear and quadratic fit of 
niBUp. 



Here we briefiy discuss the effect of one-loop correction to these quantities. Figure |3 



shows the renormalized values of V4^{iq = 1), Uk{iq = 2), and niBUpiiq = 1) at k = 0.1570. 
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As mentioned at the end of Section |T|, the leading logarithmic factor Eq. ([TTD is multiplied to 
i^. We also list the results of linear fits of them in Table |^ As we discussed previously, the 
l/mg dependence of the one-loop coefficient is significant only for V4 and almost negligible 
for Vi- As a result, the l/m^ dependence of V4 is largely affected by the renormalization 
effect, and it even changes the sign of the slope in l/m^. The I/tub dependence of V4 is still 
mild after the renormalization effect is included. For f/^ and msUp the l/rriB dependence 
is not affected by the one- loop correction, while their amphtudes decrease by at most 30%. 



B. g^-dependence of the form factors 

First we study for which region our present statistics allow us to compute the form 
factors with reasonable statistical errors. The dependence of the form factors and 
/° are shown in Figs, n and at mQ=2.6 and 1.5, respectively. We find that for k = 



0.1570(mg ~ 2ms), the range of in which the form factors have good signal covers almost 
the entire kinematic region for D meson and one third of the kinematic region for B meson. 
For K = 0.1600 (mg ~ rUs), the signal becomes much noisier, but still the form factors 
have marginally good signal for half and one fourth of the kinematic region for D meson 
and B meson, respectively. Although our present results are very noisy after the chiral 
extrapolation, this will be improved by future high statistics studies. This is encouraging in 
view of the fact that the future B Factories can produce 10*^ B-B pairs and the branching 



fraction of B ^ nlu from CLEO is (1.8 ± 0.4 ± 0.3 ± 0.2) x 10~^ [Q. It is reasonable to 
expect that there is a possibility of observing B nlu events in the regime which the 
present lattice calculation can cope with. 

Secondly we study the dependence to see whether the contribution from the B* res- 
onance to the form factor can actually be observed in the simulation data. At the chiral 
limit, unfortunately, the results are too noisy to discuss their dependence, therefore we 
use the finite mass results only in the following analysis of the g^ dependence. As shown in 



Figs. |I4| and 0, the lattice results are available only in the large g region, at which the 
recoil momentum of pion is small enough. Therefore it is justified to express the functional 
form of the form factors by an expansion around the zero recoil limit. For this purpose we 
use the inverse form factors l//'''(g^) and l//°(g^): 

l/f{q') = l/fiqlJ + ci(gL - q') + C2(gL - Q'r- (35) 



Figure |16| shows the inverse form factors at mg = 2.6 as well as their fitted functions with 
this form. The numerical results of the fit with and without the condition C2 = are given 
in Table ^ for mq = 2.6, 1.5, and 0.9. 

The pole dominance model corresponds to a special case C2 = 0, which seems to describe 
the data very well as shown in Fig. |16|. The mass of the intermediate state is given by 
'^poic — + l/(ci/^(Q'iax))5 whlch correspouds to the vector (5*) meson mass in the 
pole dominance model. Precisely speaking, the more consistent analysis is to impose the 
condition rrip^i^ = rriB* for the fit by Eq. (PB]). This constrained fit is shown with the long 
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dashed line in Fig. |T6|. It is found that now the fit do not quite agree with the data, but 
the deviation is about 10 %. 

In Fig. |1^ we also compare rrip^i^ and the measured vector meson mass as a function of 
ttlb- Again we find that there is a discrepancy between mp^ic from the unconstrained fit and 
the measured ms* , which is around few hundred MeV. Nevertheless, it is remarkable that the 
deviation remains the same order and the mass dependence of rrip^i^ has the same trend with 
niB*- We have not yet understood whether the above discrepancies can be explained from 
the remaining systematic errors such as the discretization error. But at least qualitatively 
judging from the size of the uncertainty in our calculation, our data is not inconsistent with 
the picture that there is a sizable contribution from the B* pole to the form factor /"*" near 

• 

^ max 

So far the discussion have been based on the tree level study. Let us now study how 
one-loop renormalization changes the form factors. Because the one-loop correction is dif- 



ferent for V4 and Vi, the shape of the form factors may change significantly. Figure shows 
the form factors for rriQ = 2.6 and n = 0.1570 with renormalization factors. The leading 
logarithmic factor Eq. (|T^ is not multiplied in the present case. We find that the renormal- 
ized has stronger dependence than that of at the tree level, while /° receives only a 
small change. The renormalization makes the B* pole fit even worse. In fact, the deviation 
of the constrained fit from our renormalized /"^ data is as large as 25 % near q"^^^. This is 
still within the typical size of 0{a) errors. It is very important to perform the analysis with 
larger (3. 



C. Soft pion theorem 

Applying the soft pion theorem to the B ^ n matrix element, f^{qt^) is related to the 
B meson decay constant |lT2|J13| , p!5[ 

fiQlJ = fs/U (36) 



in the massless pion limit. This relation is examined in Fig. |T9|. For the values of /b, we 
refer our work on Jb [0] , which is obtained with an evolution equation of a slightly different 
form from that of the present work. We observe a large discrepancy between and the 
decay constant both for the l/m^ dependence and for the value itself, fs increases rapidly 
toward heavier heavy quark masses, while f^iq^^J) almost stays constant. 

The discrepancy still remains significant when the renormalization effect is incorporated. 
In evaluating the renormalized values of Jb, we use one- loop perturbative coefficient obtained 
in the same manner as in Section [II C| . The leading logarithmic factor Eq. (|1^ is multiplied 
to both fiqiJ and /b. 

One may argue that the observed discrepancy can be explained by the uncertainty in the 
extrapolation procedure. To study this possibility, we compare f^{qt^J and fs/ also in 
finite light quark mass cases, in the light of the heavy meson effective theory which implies 
the relation (pSl) . They are compared in Fig. ^ as a function of 1/k. The difference between 
them are remarkable even for finite light quark mass cases. 
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The reason why these differences occur is not clear. Since our present results suffer from 
various systematic uncertainties, as described in the next section, further study with better 
control of systematic errors is necessary to clarify the origin of the problem. 



V. SYSTEMATIC ERRORS 

In this section, we qualitatively discuss on the systematic uncertainties associated with 
the lattice regularization. The following is a list of the main sources of systematic errors: 

• 0(a) errors: The characteristic size of 0{aKQCD) error arising from the unimproved 
Wilson quark action at /? = 5.8 is 20-30%. This effect is large enough to explain the 
discrepancy between EsijPjUp + ET,{k)Uk and V^, mentioned in Section |T|. Use of the 
0(a)-improved Clover action for the light quark will reduce this error to the level of 5 
%. 

• 0((ap)^) error: The systems with finite momentum may suffer from the discretization 
errors more seriously than that at the zero recoil point. The analytic estimate of the 
momentum dependent error |T^ shows that the effect is about 20 % at |p| ~ 1 GeV 



even one uses the 0(a)-improved current. 

Perturbative corrections: The one-loop correction could become significant especially 
for small (3 values. Strictly speaking, our calculation does not treat the one-loop 
effects correctly, because the stabilizing parameter n does not have correct values. 
This problem must be removed in the future studies. In estimating the one-loop 
corrections, we did not include the effect of the operator mixing, which was reported 
to be significant in the case of fs [0]- This effect also should be included to obtain 



reliable results. 

• Oil/mq) effects: We described the heavy quark with the NRQCD action including 
the order l/mq terms. Further precise calculations may need to include ©(l/mg) 
corrections, although the effect was shown to be small [0 for Jb- 

The finite volume effect may also be important. 

Since the all above systematic errors can be large, there is no advantage of giving quan- 
titative estimates of each error at this stage. The use of the 0(a)-improved (clover) action 
for light quark, as well as the simulation at higher (3 values will reduce most of the above 
systematic errors. The simulation with dynamical quarks is also of great importance for 
reliable predictions of the weak matrix elements. 

VI. CONCLUSION 

In this paper, we present the results of the study of S ^ vr form factors using NRQCD 
to describe the heavy quark with the Wilson light quark. Clear signal is observed for the 
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matrix element in a wide range of heavy quark mass containing the physical 6-quark mass. 
They are extrapolated to the chiral limit, although the result is so noisy for quantitative 
conclusion. 

The l/rriB dependence of the matrix elements are studied and it is clarified that the 
temporal component and the part of the spatial component proportional to the pion mo- 
mentum have fairly small dependencies on mq. On the other hand, the part of the spatial 
component proportional to the B momentum has a significant ©(l/ms) correction. 

The dependence of the form factors in the finite light quark masses are studied. We 
find that the dependence of the form factor f^{q'^) near q^^^^ becomes much stronger for 
larger heavy quark mass. Model independent fit of 1/ f^{q^) near q^^^^ shows that the tree 
level results are consistent with the pole behavior for large q^ range, and the difference of 
fitted pole mass and the measured m^* is around few hundred MeV for all the heavy quark 
masses. 

The values of /° at the zero recoil point are compared with the prediction of the soft 
pion theorem, and the significant discrepancy is observed. 

The size of the renormalization corrections are estimated by the one-loop perturbative 
calculation. They almost does not affect their l/mq dependence, but decrease Vi much 
more than V4, which drastically change the shape of Our present result suffers from 
large systematic uncertainties, and the most important one is 0{a) error. It is very important 
to study at higher (5 with improved actions. 
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TABLES 

(mQ.n) A D Cy, Cy 

(5.0,1) 0.0759 0.0124(4) 0.0210(11) -0.0790(10) 

(2.6,2) 0.0668 0.0353(3) 0.0004(9) -0.0780(7) 

(2.1.2) 0.0623 0.0449(3) -0.0068(9) -0.0757(7) 

(1.5.3) 0.0528 0.0623(2) -0.0192(8) -0.0734(6) 
(1.2,3) 0.0446 0.0757(1) -0.0283(8) -0.0707(6) 
(0.9,6) 0.0309 0.0933(1) -0.0428(8) -0.0687(5) 



TABLE I. The tadpole improved one- loop coefficients for the perturbative corrections Eq, 
Zjn, Zy^, and Zy^- Quoted errors represent the numerical uncertainties in the evaluation of loop 
integrals. The uncertainty of A is less than 10~^. 



At = 0.1570 0.1585 0.1600 Kc 

0.5677(30) 0.4933(33) 0.4118(37) 

nip 0.6747(54) 0.6214(72) 0.567(11) 0.448 (17) 

U 0.1496(46) 0.1380(49) 0.1270(53) 0.1019(64) 

TABLE IL The values of m,r, rnp, and pion decay constant without renormalization. Fitting 
range is t = 14 — 24. 
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TABLE III. The momentum combinations {p, k) used in the simulation. In this table, the 
values of p, k, and q are expressed in the unit of 27r/16. The set which is equivalent with another 
under the lattice rotational symmetry is identified with the same ig-number, and a representative 
is shown in the fifth through seventh columns. The last column shows the numbers of equivalent 
combinations. The symbols in the third column denote the direction of k against p as follows: 
±: orthogonal, ft- parallel, jj: anti-parallel, and oblique for others. The set ig = 12 gives the 
minimum value among the sets in this table. 
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Pseudosealar meson l)inding energy: E,i(j(p = 0) 



("^Q , n) 


K = 0.1570 


0.1585 


0.1600 




(5.0,1) 


0.6304(69) 


0.6084(83) 


0.585 (11) 


0.535 (15) 


(2.6,1) 


0.6268(48) 


0.6041(56) 


0.5809(71) 


0.530 (10) 


(2.1,1) 


0.6247(45) 


0.6014(52) 


0.5777(65) 


0.5260(91) 


(2.1,2) 


0.6279(53) 


0.6056(62) 


0.5834(80) 


0.534 (11) 


(1.5,2) 


0.6180(42) 


0.5940(48) 


0.5696(59) 


0.5162(81) 


[L.Z ,Z) 


U.DioO(^4U j 


U.0ooy(4D j 


U.0D4U(0D j 


U.OUyo( iO) 


(1.2,3) 


0.6142(51) 


0.5899(56) 


0.5655(68) 


0.5117(92) 


(0.9,2) 


0.6058(39) 


0.5805(43) 


0.5551(51) 


0.4991(68) 






iiicoUii uiiiU-iii^ tJiici^y. 


^qQ* \P — 




{mq , n) 


K = 0.1570 


0.1585 


0.1600 




(5.0,1) 


0.649 (12) 


0.628 (14) 


0.604 (19) 


0.555 (27) 


(2.6,1) 


0.6502 (62) 


0.6287 (76) 


0.6065 (99) 


0.559 (14) 


C9 1 1 ^ 


U.OOUl ^^OO ) 








(1.5,2) 


0.6488 (52) 


0.6257 (61) 


0.6014 (79) 


0.550 (11) 


(1.2,2) 


0.6484 (51) 


0.6249 (59) 


0.6002 (76) 


0.547 (11) 


(0.9,2) 


0.6470 (50) 


0.6231 (57) 


0.5982 (73) 


0.545 (10) 



TABLE IV. The binding energy of the pseudosealar and vector heavy-hght mesons. The single 
exponential fit is applied with the fitting range t = 1(S — 24. For the pseudosealar we average the 
values obtained from the local-local and the smeared-local correlation functions. For the vector 
mesons we use the local-local only, and there is no data available for (jtzq, n) = (2.1, 2) and (1.2, 3). 
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<f 


% 


Up 


Uk 


1 


7.071 (20) 


1.014 (34) 


- 


- 


2 


6.280 (19) 


0.844 (26) 


- 


0.878 (41) 


3 


5.609 (19) 


0.754 (50) 


- 


0.695 (61) 


4 


5.017 (18) 


0.612 (87) 


- 


0.57 (10) 


5 


7.044 (20) 


0.999 (36) 


0.0475(28) 


- 


6 


6.247 (19) 


0.832 (28) 


0.0366(47) 


0.860 (41) 


7 


6.555 (19) 


0.930 (30) 


1.009 (46) 


1.009 (46) 


8 


5.938 (19) 


0.750 (34) 


-0.702 (48) 


0.702 (48) 


9 


5.571 (19) 


0.742 (49) 


0.040 (12) 


0.674 (59) 


10 


5.880 (19) 


0.827 (55) 


0.790 (68) 


0.767 (66) 


11 


5.283 (18) 


0.66 (10) 


0.65 (12) 


0.63 (11) 


12 


4.666 (18) 


0.544 (68) 


-0.39 (12) 


0.477 (82) 


13 


7.017 (20) 


0.992 (42) 


0.0467(30) 


- 


14 


6.214 (19) 


0.825 (34) 


0.0360(48) 


0.848 (45) 


15 


6.523 (19) 


0.923 (38) 


0.517 (26) 


0.997 (51) 


16 


5.534 (19) 


0.757 (76) 


0.052 (53) 


0.670 (82) 


17 


6.151 (19) 


0.920 (67) 


0.863 (77) 


0.863 (77) 


18 


5.842 (19) 


0.820 (58) 


0.412 (36) 


0.758 (68) 


19 


5.225 (19) 


0.669 (52) 


-0.266 (41) 


0.587 (61) 


20 


6.990 (20) 


0.968 (58) 


0.0454 (33) 





TABLE V. V4, Up, and Uk in the lattice unit at mg = 2.6 and k = 0.1570. iq denotes the 
set of momentum {p, k) summarized in Table III , In the evaluation of , the B meson mass is 
determined through the tree level relation tub = rnq + EqQ(0). 
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linear 




quadratic 
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4 


4 


4 


4 


4 




1 




1 RA(Fi?\) 










2 


0.826(29) 


0.080(47) 


0.851(41) 


-0.06(17) 


0.15(17) 




3 


0.757(51) 


-0.038(59) 


0.799(57) 


-0.30(20) 


0.31(22) 




4 


0.624(80) 


-0.25 (11) 


0.79 (10) 


-1.29(36) 


1.25(42) 




1 




165 ('631 


1 ni6f55'l 


—0 0nf23'l 






2 


0.807(35) 


0.075(57) 


0.830(48) 


-0.06(20) 


0.14(19) 




3 


0.758(76) 


-0.071(73) 


0.830(81) 


-0.51(26) 


0.51(29) 




4 


0.62 (12) 


-0.40 (15) 

V / 


0.89 (19) 


-1.83(50) 


1.75(60) 


0.1600 


1 


1.003(53) 


0.150(76) 


1.023(66) 


0.05(27) 


0.10(25) 




2 


0.768(46) 


0.088(76) 


0.788(58) 


-0.04(26) 


0.14(25) 




3 


0.78 (14) 


-0.17 (10) 


0.96 (17) 


-1.13(40) 


1.13(46) 




4 


0.70 (27) 


-0.64 (25) 


1.22 (55) 


-2.45(80) 


2.26(94) 



TABLE VI. Parameters for the linear and quadratic fits of V4{p = 0,k). 
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linear quadratic 



K 


ig 










^k 


0.1570 


2 


0.945(39) 


-0.194(44) 


0.967(47) 


-0.30(19) 


0.13(19) 




3 


0.762(56) 


-0.257(53) 


0.750(54) 


-0.17(22) 


-0.10(24) 




4 


0.655(88) 


-0.364(91) 


0.600(81) 


0.08(43) 


-0.54(49) 


0.1585 


2 


1.004(52) 


-0.198(50) 


1.023(58) 


-0.28(22) 


0.10(23) 




3 


0.808(92) 


-0.242(64) 


0.769(80) 


0.00(30) 


-0.29(32) 




4 


0.72 (15) 


-0.34 (14) 


0.58 (12) 


0.77(74) 


-1.34(80) 


0.1600 


2 


1.064(73) 


-0.214(62) 


1.063(77) 


-0.21(29) 


0.00(30) 




3 


0.92 (20) 


-0.219(90) 


0.80 (16) 


0.47(50) 


-0.81(53) 




4 


0.94 (37) 


-0.23 (26) 


0.55 (23) 


3.3 (2.3) 


—A 1 O 




TABLE VII. Parameters for the linear and quadratic fits of Uk{p = 


U, K). 








linear 




quadratic 




K 


iq 








Cp 


„(2) 


0.1570 


1 


0.0887(80) 


2.61(39) 


0.0717(95) 


4.5(1.2) 


-1.55(76) 




2 


0.089 (14) 


1.29(38) 


0.072 (13) 


2.7(1.1) 


— 1.31(8oj 


0.1585 


1 


0.0872(94) 


2.65(47) 


0.066 (11) 


5.3(1.7) 


-2.1 (1.0) 




2 


0.093 (20) 


0.98(42) 


0.080 (18) 


2.0(1.2) 


-1.0 (1.1) 


0.1600 


1 


0.088 (12) 


2.72(59) 


0.059 (15) 


6.7(2.7) 


-3.1 (1.7) 




2 


0.104 (33) 


0.67(47) 


0.097 (27) 


1.1(1.5) 


-0.4 (1.5) 



TABLE VIII. Parameters for the linear and quadratic fits of Up{p = 0,k). 
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0.1570 
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0.052(55) 


1.088(39) 




-0.209(47) 


0.1585 


1.019(44) 






0.039(63) 


1.105(46) 




-0.216(55) 


0.1600 


1.0.39(55) 






0.030(77) 


1.126(58) 




-0.219(66) 








UkKP ~- 


n 1 L 1 1 ^ 
— \J^\K\ — L) 


( «g = 2 ) 








MO) 


q* 


= vr/a 


(1) 




q* = 1/a 


(1) 


0.1570 


0.732(31) 






0.013(61) 


0.609(27) 




0.081(70) 


0.1585 


0.778(42) 






0.005(68) 


0.649(36) 




0.070(78) 


0.1600 


0.826(59) 






-0.019(84) 


0.689(50) 




0.043(96) 










— 0, K — 0) 


I - 1 J 






K 


7-7(0) 
Up 


q* 


= 7r/a 


Cp 


Up 


q* = 1/a 




0.1570 


0.0466(66) 






6.3(1.2) 


0.0268(58) 




11.3(30) 


0.1585 


0.0453(77) 






6.5(1.5) 


0.0256(68) 




11.8(38) 


0.1600 


0.045 (10) 






6.7(1.9) 


0.0248(87) 




12.5(53) 



TABLE IX. Parameters for the linear fit of the renormalized matrix elements V4^{p = 0,k = 0), 
Uk{p = 0,\k\ = 1), and Up{p = 0,k = 0). 
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hnear fit 




quadratic fit 




{mQ,n) K 




Cl 




Cl 


C2 


(2.6, 1) 0.1570 


1.373(54) 


0.126(70) 


1.386(52) 


0.058(64) 


0.046(53) 




0.480(21) 


0.264(38) 


0.470(20) 


0.335(40) 


-0.051(37) 


0.1585 


1.436(70) 


0.109(88) 


1.438(64) 


0.098(91) 


0.007(81) 




0.445(24) 


0.272(47) 


0.434(22) 


0.366(59) 


-0.068(54) 


0.1600 


1.531(94) 


0.09 (11) 


1.512(86) 


0.22 (16) 


-0.09 (14) 




0.407(27) 


0.276(61) 


0.395(26) 


0.44 (10) 


-0.115(86) 


(1.5,2) 0.1570 


1.167(38) 


0.209(81) 


1.185(37) 


0.086(80) 


0.119(87) 




0.597(25) 


0.472(64) 


0.587(22) 


0.548(60) 


-0.075(78) 


0.1585 


1.213(50) 


0.19 (10) 


1.224(47) 


0.10 (12) 


0.08 (14) 




0.559(28) 


0.493(78) 


0.545(24) 


0.623(92) 


-0.13 (12) 


0.1600 


1.283(67) 


0.17 (14) 


1.279(62) 


0.21 (20) 


-0.04 (24) 




0.516(32) 


0.52 (10) 


0.496(29) 


0.77 (17) 


-0.26 (20) 


(0.9,2) 0.1570 


1.011(28) 


0.360(88) 


1.027(27) 


0.208(85) 


0.19 (13) 




0.685(28) 


0.753(90) 


0.690(26) 


0.713(75) 


0.05 (13) 


0.1585 


1.041(36) 


0.35 (11) 


1.056(35) 


0.19 (13) 


0.20 (21) 




0.647(33) 


0.79 (11) 


0.640(28) 


0.86 (12) 


-0.09 (21) 


0.1600 


1.090(49) 


0.33 (15) 


1.096(48) 


0.26 (24) 


0.10 (35) 




0.599(37) 


0.85 (14) 


0.577(32) 


1.12 (23) 


-0.36 (36) 


TABLE X. Parameters for the fit f-^{q^) = /" 




-g2) + c2(gi^ 


— q^Y-, where 


C2 is set to zero for the hnear fit. 


For each (mQ, 


n) and n, numbers in upper and lower rows 


correspond and 


respectively. In ah the cases 


, x^/dof are less than unity. 
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FIG. 1. Renormalization constant for the vector current with two scales for the couphnj 
constant, q* = ir/a and 1/a. The open and filled symbols represent Zy^ and Zy., respectively. 
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FIG. 2. Effective mass plot of pion at k = 0.1570 and 0.1600. The horizontal sohd hnes 
represent the fitted values and the fitting range with the statistical errors (dotted lines). 
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FIG. 3. Dispersion relation for pion. The solid lines represent the relation E^{k) = + /c^ 
with the rest mass obtained in the simulation. For k = 0.1585 and 0.1600, symbols are slightly 
shifted in horizontal direction for clarity. 
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FIG. 4. Effective mass plot of B meson at niQ = 2.6 and k = 0.1570, 0.1600. Results with the 
smeared source (filled symbols) are shown for = 0, 1 as well as the results with the local source 
(open symbols). The horizontal solid lines express the average values over the results of single 
exponential fit of the local-local and the smeared-local correlation functions. The statistical errors 
of the fitted values are displayed at the right end of the lines. For all niQ, k, and momentum, the 
fit ranges are set to t = 16 — 24. 
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FIG. 5. Dispersion relation for the B meson at mq = 2.6 and k = 0.1570, 0.1600. The sohd 
hnes represent the relation EqQ{p) = EqQ{0) + /2mB, for which is determined with the tree 
level formula ruB = mq + EqqiG). Dashed lines represent the same relation with the renormalized 
at the scale q* = 1/a. 
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FIG. 6. Effective mass plot for the three-point functions at mq = 2.6 and k = 0.1570, 0.1600. 
The horizontal lines express the values obtained from the two-point correlation functions with the 

(3) 

statistical errors indicated at the right end of the lines. Top two figures are for C4 and the third 
is for fc-C(^) at AC = 0.1570, and bottom figure is for cf ^ at k = 0.1600. 
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FIG. 7. Comparison of % to ( EsiTp + E^Uk ) for = 6 at k = 0.1570. 
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FIG. 8. Chiral extrapolation of the matrix elements for mq = 2.6. V4 and k-V/\k\^ are shown 
for two momentum configurations iq = 1, 2. The solid and the dashed lines represent the linear fit. 
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FIG. 9. V4 at K = 0.1570. For k = and \k\ = 1, the soHd and the dashed hnes represent the 
results of linear and quadratic fits, respectively. 
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FIG. 10. Uk at K = 0.1570. The solid and the dashed lines represent the results of linear and 
quadratic fits, respectively. 
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FIG. 11. Extraction of Up{p = 0,k) is shown for uiq = 2.6 and 1.5 at k = 0.1570. The 
extrapolation is carried out hnearly in p^. For k = 0, iq = 5, 13, 20 are used. For k = 1, Up{0, k) is 
determined using iq = 6, 14, for which p and k are perpendicular. 
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FIG. 12. Up multiplied by at k = 0.1570 as a function of l/ms. The values of me are 
determined with the tree level formula. The solid and the dashed lines represent the linear and 
the quadratic fits, respectively. 
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FIG. 13. One-loop renormalized V4, Uk-, and tip as a function of l/mg. The solid, the dashed, 
and the long dashed lines represent the linear fits. 
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FIG. 14. Form factors at mg = 2.6 and k = 0.1570, 0.1600. The solid curves represent the fit 
to single pole functions. 
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FIG. 15. Form factors at mg = 1.5 and k = 0.1570, 0.1600. 
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FIG. 16. l/f and 1//+ as a function of at mq = 2.6 and k = 0.1570. The solid and 
the dotted curves represent the linear and the quadratic fits, respectively. The long dashed line 
represents the linear fit with the constraint nZpoie = itib*, where m*^ is the B* meson mass obtained 
from the two-point correlation function. 
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FIG. 18. Renormalized form factors at mq = 2.6 and k = 0.1570. Upper and lower figures 
are obtained with gy{7r/a) and gy{l/a), respectively. The solid lines represent the results of the 
linear fit. For f~^{q'^), x^/dof are 1.9 and 4.1 for q* = tt/o and 1/a respectively, x^/dof for f^{q'^) 
are less than 0.5 in the both cases. 
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